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ABSTRACT 

We use deep near-infrared and submillimeter observations of three massive lensing cluster fields, A370, 
A851, and A2390, to determine the average submillimeter properties of a i^'-selected sample. The 38 
Extremely Red Objects (EROs; I — K' > 4) with K' < 21.25 have a very significant error-weighted 
mean 850 /j,m flux of 1.58 ± 0.13 mJy. The ERO contribution to the 850 /im background is 1.88 ± 0.16 x 
10 4 mJy deg" 2 , or about half the background light. The 17 Very Red Objects (VROs; 3.5 < I- K' < 4) 
are also significantly detected (1.32 ± 0.19 mJy), bringing the combined VRO and ERO contribution to 
2.59 ± 0.19 x 10 4 mJy deg" 2 . There is a substantial systematic uncertainty (about a factor of two) in 
this value due to field-to-field variation, but even with this uncertainty it is clear that a large fraction 
of the 850 fim background arises from red objects. An analysis of the VRO and ERO number counts 
shows that half of the population's 850 light arises in objects with demagnified magnitudes K' < 20 
and half in fainter objects. On the basis of the I — J versus J — K' colors of the galaxies, the bulk of 
the submillimeter signal appears to arise from the dusty starburst galaxies in the red object population 
rather than from the high-redshift elliptical galaxies. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Extremely Red Objects (EROs) have been the subject 
of many studies since their discovery (e.g., Cowie et al. 
1990; Hu & Ridgway 1994). Two scenarios can explain 
the unusally red colors of the EROs: old stellar popula- 
tions at z > 1 (e.g., Dunlop et al. 1996; Spinrad et al. 1997; 
Soifcr et al. 1999; Cimatti et al. 2002) or dust-enshrouded 
galaxies (e.g., Graham & Dey 1996; Cimatti et al. 1998; 
Smith et al. 2001). The starlight from dusty galaxies is re- 
processed and reemitted in the rest wavelength far-infrared 
band. For high redshift sources, this light is redshifted into 
the submillimeter. Recently, a number of possible ERO 
counterparts to submillimeter detected sources have been 
discovered (e.g., Smail et al. 1999, 2002; Barger, Cowie, & 
Richards 2000; Ivison et al. 2002). 

The COBE satellite found that the submillimeter extra- 
galactic background light (EBL) has approximately the 
same integrated energy density as the optical EBL (Puget 
et al. 1996; Fixsen et al. 1998); thus, a large fraction of 
starlight is reradiated by dust. To determine the high- 
redshift star formation, the sources that comprise the sub- 
millimeter EBL need to be determined. Here we use deep 
submillimeter and near-infrared (NIR) imaging of three 
massive lensing cluster fields to determine the submillime- 
ter properties of the ERO population. Two advantages 
of observing cluster fields are that background sources are 
magnified and mean source separations on the sky are in- 
creased. 

In § |^ we briefly describe our /, J, K', and 850 /j,m 
data. The 15" beam size of SCUBA (Holland et al. 1999) 



on the 15 m James Clerk Maxwell Telescope 5 prohibits 
us from making individual matches between submillime- 
ter and optical/NIR sources. However, in § || we analyze 
our data statistically and find that EROs, as a class, have 
a significant submillimeter flux and produce much of the 
850 fim EBL. In § [I] we use our NIR photometry and lens- 
ing corrections to calculate the ERO surface densities to 
K' = 24. We analyze where in K 1 magnitude the 850 /im 
flux arises. Using the /, J, K' colors, we differentiate 
between evolved galaxies and dust-reddened starbursts to 
show that the 850 /jm light primarily arises in the latter 
class of object. 



2. DATA AND SAMPLE SELECTION 

We constructed K' catalogs for our three cluster fields 
using data obtained with the Cooled Infrared Spectro- 
graph and Camera for OHS (CISCO; Motohara et al. 
1998) on the Subaru 6 8.2 m telescope. The detector is 
a 1024 x 1024 HgCdTe HAWAII array with a 0.111" pixel 
scale and an ~ 2' x 2' field-of-view. The K' exposures 
are 15.4 ks (A370), 3.84 ks (A851), and 15.4 ks (A2390). 
J-band data were also obtained with the same instru- 
ment. Details of the observations and data reduction can 
be found in Cowie et al. (2001; A370 and A2390) and 
L. L. Cowie et al., in preparation (A851). The image 
quality of the K 1 images is 0.66" FWHM. We used K' 
magnitudes measured in 2" diameter apertures to deter- 
mine colors, and we used isophotal magnitudes to 2% of 
the peak surface brightness for total magnitudes. We es- 

1 Dcpartment of Astronomy, University of Wisconsin-Madison, 475 North Charter Street, Madison, WI 53706 
2 Department of Physics and Astronomy, University of Hawaii, 2505 Correa Road, Honolulu, HI 96822 
^Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, Hawaii 96822 
4 Obscrvatoire Midi-Pyrenees, 14 Avenue E. Belin, 31400, Toulouse, France 

5 The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre on behalf of the UK Particle Physics and Astronomy 
Research Council, the Netherlands Organization for Scientific Research, and the Canadian National Research Council 
6 The Subaru telescope is operated by the National Astronomical Observatory of Japan. 
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timated the noise in the images by measuring aperture 
magnitudes at a number of blank sky positions. The 5a 
limit for the shallowest of the fields is K' — 21.25; we 
adopt this limit for the entire K' sample selection. 

We measured 2" diameter aperture /-band magnitudes 
for the K' samples from images obtained with either the 
Low- Resolution Imaging Spectrometer (LRIS; Oke et al. 
1995) or the Echellette Spectrograph and Imager (ESI; 
Epps & Miller 1998) on the Keck 7 10 m telescopes. The 
data are described in Cowie et al. (2001; A370 and A2390) 
and Smail et al. (1999; A851). The 2a limit for the shal- 
lowest of the fields is / = 26; we adopt this limit for the 
entire sample. 

Finally, we measured 850 /jm fluxes at the K' source 
positions from ultradeep SCUBA jiggle maps. We used 
beam-weighted extraction routines that include both the 
positive and negative portions of the beam profile. The 
total exposure times are 133.1 ks (A370), 59.3 ks (A851), 
and 78.1 ks (A2390). The la sensitivities, including both 
fainter sources and correlated noise, are 0.45 mJy (A370), 
0.87 mJy (A851), and 0.71 mJy (A2390). The data are 
described in Cowie, Barger, & Kneib (2002). 

3. SUBMILLIMETER PROPERTIES OF EROS AND THEIR 
CONTRIBUTION TO THE SUBMILLIMETER EBL 



high-redshifL evolved cluster lying behind A2390, as we 
discuss in § Q. 
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Fig. 2. — (a) 850 (im flux and (b) error-weighted mean 850 fim 
flux versus J — K' color for the combined K' samples in the A370, 
A851, and A2390 fields. In (a) any sources with 850 /im flux greater 
than 9.5 mJy are plotted at 9.5 mJy. In (b) the vertical bars show 
the 1<t uncertainties and the horizontal bars show the color range 
sampled in each bin. 

The 15" SCUBA beam size means that we cannot be 
certain which if'-selected galaxy is the true counterpart 
to a submillimeter source. However, by determining the 
error-weighted mean 850 /im flux of the K' sources, we 
can make statistically significant statements regarding the 
relationship between submillimeter flux and source color. 
This is analogous to what Peacock et al. (2000) did for 
optical starbursts. 

In Fig. ||a we plot submillimeter flux versus I — K' color, 
and in Fig. IJd we plot error-weighted mean 850 /im flux per 
color bin. Because submillimeter sources are not removed 
from the jiggle maps, the complex beam profile produced 
by the nod and chop observing pattern (there are twice as 
many negative as positive response points, but at half the 
sensitivity) means that the submillimeter flux distribution 
is non-Gaussian with an extended negative tail. A nega- 
tive submillimeter flux will be measured for a iC'-selected 
source positioned in the negative portion of the beam pro- 
file corresponding to a submillimeter source. Thus, there 

7 The W. M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of 
California, and NASA, and was made possible by the generous financial support of the W. M. Keck Foundation. 
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Fig. 1. — I — K' versus K' for the combined K' samples in 
the A370, A851, and A2390 fields. The vertical line denotes the 
K 1 = 21.25 (5<r) limit adopted. The upward pointing arrows denote 
sources that were not detected above I = 26 (2a). 

Figure [l] shows the I — K' versus K' color-magnitude 
diagram. Most of the sources are cluster members. The 
cluster color-magnitude relationship is clearly visible and 
matches the measured relations for clusters at these red- 
shifts (e.g., Kodama et al. 1998). Thirty-eight galax- 
ies have I — K' > 4, classifying them as EROs, and 
fifty-five have I — K' > 3.5. We classify galaxies with 
3.5 < / - K' < 4 as Very Red Objects (VROs). The 
11.5 arcmin 2 area covered by the K' data is somewhat 
smaller than the 18.9 arcmin 2 area covered by the submil- 
limeter data. A disproportionate number of the EROs 
and VROs lie in the A2390 cluster field (26 of the 38 
sources with I — K' > 4, and 32 of the 55 sources with 
I — K' > 3.5). The excess may be associated with a very 



Wehner, Barger, & Kneib 



3 



are more high negative flux measurements in Fig. ||a than 
would be expected based on a simple Gaussian model with 
the quoted errors. The average signal of a set of random 
sources should, however, be zero, and the significance of 
the result can be determined from sets of Monte Carlo 
simulations of the same number of points randomly dis- 
tributed in the fields (more details can be found in Cowie 
et al. 2002). 

At low values of J — K' the error- weighted mean 850 /mi 
fluxes are consistent with zero; however, as I — K' in- 
creases, there is a corresponding increase in the means 
(see Fig. |^b). The I — K' > 4 galaxies have an error- 
weighted mean of 1.58 ± 0.13 mJy. The formal error is 
based on the assumption of Gaussian statistics and may 
underestimate the true error, but the result is highly sig- 
nificant: 10000 Monte Carlo simulations yielded only one 
case where the same number of randomly positioned ob- 
jects had a value equal or above the observed value. We 
tested whether the signal is arising from a small number 
of objects by eliminating the positive and negative tails: 
if the two (five) most positive and two (five) most nega- 
tive sources are removed from the / — K' > 4 sample, the 
signal becomes 1.40 ± 0.14 mJy (1.48 ± 0.15 mJy). The 
I — K' > 4 objects in the A2390 field have a value of 
1.41 ± 0.18 mJy, while those in the A370 and A851 fields 
have a value of 1.80 ±0.20 mJy, so the mean submillimeter 
flux per ERO is consistent, even if we break down the sam- 
ple. The 3.5 < I — K' < 4 galaxies have an error-weighted 
mean of 1.32 ± 0.19 mJy; here 10000 simulations yielded 
only 85 cases with values at or above this level. 

Since the sky surface brightness is conserved by the 
lensing process, we find the ERO population contributes 
1.88 ± 0.16 x 10 4 mJy deg" 2 to the 850 /on EBL. This 
value is 61 percent of the total 850 /mi EBL if the 
Puget et al. (1996) value of 3.1 x 10 4 mJy deg -2 is as- 
sumed, or 43 percent if the Fixsen et al. (1998) value of 

4.4 x 10 4 mJy deg -2 is assumed. When the galaxies with 

3.5 < / — K ' < 4 are included, the contribution increases 
to 2.59 ±0.19 x 10 4 mJy deg~ 2 (83 or 59 percent of the to- 
tal, respectively). Thus, the majority of the 850 /im EBL 
arises from the relatively small number of sources with 
I - K' > 3.5. 

However, there may be systematic errors arising from 
clustering. For the individual fields, the EBL contribu- 
tions from the I — K' > 3.5 sources are 0.81 ±0.14 x 
10 4 mJy deg" 2 (A370), 6.82±0.42xl0 4 mJy deg~ 2 (A851), 
and3.49±0.55xl0 4 mJydeg- 2 (A2390). Although A2390 
has the largest number of / — K' > 3.5 objects (see § Q), 
A851 has a higher mean submillimeter signal per object. If 
we exclude A2390 because of the anomalously high number 
of red objects in this field, then the I — K' > 3.5 popula- 
tion contributes 1.94±0.15 x 10 4 mJy deg -2 to the 850 /im 
EBL. While the I — K' > 3.5 populations in all three clus- 
ter fields are significantly detected, it is clear that a large 
number of fields would have to be averaged to produce a 
high precision estimate of the EBL contribution, and the 
systematic errors in the present estimate may be as high 
as a factor of 2. 

4. THE BREAKDOWN OF THE ERO CONTRIBUTIONS 

We now analyze the ERO number counts to determine 
how representative our fields are and to estimate the ERO 
surface densities to fainter magnitudes than has previously 



been possible. The methodology used is described in more 
detail in Smith et al. (2002), who analyzed 10 cluster fields 
observed to shallower K limits. 

We used the LENSTOOL software package to determine 
the amplifications of the EROs, assuming they were lo- 
cated at z = 2. (The amplification varies only weakly with 
rcdshift for any source beyond z = 1 that has a modest am- 
plification.) LENSTOOL uses multiple-component mass 
distributions that describe the extended potential wells of 
the clusters and their more massive individual member 
galaxies (Kneib et al. 1996). Details for the A370, A851, 
and A2390 clusters can be found in Kneib et al. (1993), 
Scitz et al. (1996), and J. -P. Kneib et al., in preparation, 
respectively. 

We also used LENSTOOL to determine the areas in the 
source plane over which an object with a given source plane 
magnitude would be seen. We created grids the size of our 
cluster field K' images with 1" x 1" grid elements and 
treated the grid points as objects in the source plane at 
z = 2. We brought the grids through the cluster mass dis- 
tributions to determine the magnifications at each point. 
Since each point represents one square arcsecond, the asso- 
ciated source plane area for a given magnitude is the sum 
of the points where an object of that magnitude could be 
detected. 

To compute the ERO surface densities, we summed the 
reciprocals of the associated areas for all the sources in 
each magnitude bin. (Multiply imaged sources are cor- 
rectly handled since the source plane area is proportionally 
increased relative to the number of images.) The results 
are shown in Fig. ^, where surface density per square de- 
gree per magnitude is plotted versus lensing-corrected (i.e 
source plane) isophotal K 1 magnitude. The fact that Fig 
flattens at faint magnitudes indicates that by K' = 24 we 
have already seen most of the K' light density in the ERO 
population. 




17 18 19 20 21 22 23 
SOURCE PLANE K' MAGNITUDE 



24 



Fig. 3. — Surface density per square degree per magnitude ver- 
sus lensing-corrected isophotal K' magnitude for the combined K' 
samples in the A370, A851, and A2390 fields with / - K' > 4 (filled 
squares). The uncertainties are Ict, based on a Poissonian distri- 
bution. The open diamonds show the surface densities when the 
A2390 field is excluded. 

At brighter K' magnitudes our surface densities are 
higher than those determined by Smith et al. (2002) and 
by previous blank field surveys (e.g., Barger et al. 1999; 
Daddi et al. 2000). To K' = 21.6 we find a cumulative sur- 
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face density of 8.0 (6.6, 9.4) EROs per square arcminute, 
where the parentheses show the la limits. This is slightly 
more than three times the Smith et al. value of 2.5 ± 0.4 
arcmin~ 2 . (Our I — K 1 = 4 criterion is very similar to their 
R — K > 5.3 criterion.) The reason for the excess lies in 
the large population of EROs in the A2390 field, where the 
surface density is 19.7 (15.6, 23.8) arcmin -2 to K' = 21.6. 
If we restrict our analysis to the A370 and A851 cluster 
fields (the open diamonds in Fig. ||), the surface density 
drops to 3.4 (2.3, 4.8) arcmin" 2 to K' = 21.6, which is 
comparable to the Smith et al. value. 

Inspection of the K' image of A2390 suggests that the 
ERO excess in this field corresponds to a cluster of galaxies 
centered on a bright red galaxy in the northwest corner of 
the image; this galaxy is an X-ray point source (Fabian et 
al. 2000) with a spectroscopic redshift of z = 1.467 (Cowie 
et al. 2001). A second X-ray source in the southwest corner 
has a similar probable spectroscopic redshift (Cowie et al.). 
If we assume that the excess of EROs is associated with 
this cluster, then the cluster contains roughly 20 EROs to 
K' = 21.6, making it an extremely large concentration of 
evolved galaxies at this redshift. 

The error-weighted mean 850 /xm flux per source de- 
creases (though slower than linearly) with increasing 
source plane K' magnitude: at mean K' magnitudes 
of (16.9, 18.8, 20.6, 22.4) we find error-weighted mean 
850 /im fluxes of (5.18±0.38, 1.55±0.29, 1.22±0.14, 0.47± 
0.20) mjy. Although submillimctcr flux depends only 
weakly on redshift, such that higher redshift sources fade 
in K' but not in 850 /im flux, sources that are fainter in K' 
luminosity may be intrinsically fainter in 850 /im luminos- 
ity. The dependence of 850 ^m flux on ERO K' magnitude 
means that brighter submillimeter surveys are more likely 
to identify ERO counterparts. 

The observed decrease in 850 /im flux with increasing 
K' magnitude, combined with the flatness of the ERO 
counts at the fainter magnitudes seen in Fig. [jj, implies 
that the largest contribution to the 850 /im EBL arises 
near the point at which the ERO number counts begin to 
flatten — around K' = 20. In the present data half of the 
850 /im light arises in sources with source plane magni- 
tudes K' < 20, and half arises in fainter sources, placing 
the median 850 /im source at K' = 20. If A2390, with 
its larger number of bright EROs, is excluded, the split 
occurs at K' — 20.3. When sources fainter than K' = 24 
are included, the split will occur at a fainter magnitude, 
but the effect will be small. 

We may also try to determine if the spectral energy dis- 
tributions (SEDs) of the sources correlate with the 850 /im 
fluxes. If the 850 /im signal arises primarily from dust- 
reddened starbursts rather than from old elliptical galax- 
ies, then we would expect objects with a smoother NIR 
SED to be the primary source of the signal. (The red 
/ — K' colors of old evolved galaxies are in large part due 
to the 4000 A break, which produces a discontinuity in 
the NIR SEDs.) As illustrated in Fig. [|, we have used the 
criterion I — J = 1.8( J— K') — 1.2 to separate redder I— J 
galaxies that follow the unevolved elliptical galaxy track 
from bluer I — J galaxies that lie in the dust-reddened 
starburst region (see Pozzetti & Mannucci 2000 for other 
color selection criteria). Sources in the starburst region 
could also be more recently formed galaxies or galaxies 
with some residual star formation. 



4 L I 1 1 1 1 1 I 1 1 ' I 




O F : 

12 3 4 

J-K' 



Fig. 4. — I — J versus J — K' for the combined K' samples in the 
A370, A851, and A2390 fields with measured J magnitudes. Sources 
with / — K' > 3.5 (/ — K' < 3.5) are denoted by large (small) solid 
squares; la uncertainties are given for the I — K' > 3.5 sources. 
The dashed, solid, and dotted curves show, respectively, the tracks 
for an unevolved Sb, Sa, and elliptical galaxy versus redshift from 
z = to z = 3. Redshifts z = 1, 2, and 3 are marked and labeled on 
the Sa track. The arrow shows the direction that reddening would 
move a galaxy using a simple model where extinction is inversely 
proportional to wavelength. The length of the arrow corresponds 
to Ay = 10. The dotted line I - J = 1.8(J - K') - 1.2 roughly 
separates objects that follow the unevolved elliptical galaxy track 
(left of the dotted line) from objects located where reddening could 
move a starburst galaxy (right of the dotted line). 

Of the 55 galaxies with I — K 1 > 3.5 in our sample, 39 
have measured J magnitudes. The error-weighted mean 
850 /im flux for all 39 objects is 0.90 ± 0.12 mjy. The 
17 (44 percent) that satisfy the elliptical galaxy criterion 
have 0.32 ± 0.18 mjy. In contrast, the 22 that lie on the 
dusty starburst side have 1.38 ±0.16 mjy. Thus, the bulk 
of the signal does appear to arise in the smaller number of 
galaxies that could be dusty starbursts. 

5. SUMMARY 

From an optical, NIR, and submillimeter survey of three 
massive lensing cluster fields we find that VROs and 
EROs mark, on average, strong submillimeter emitters. 
The I — K' > 3.5 population contributes 1.94 ± 0.15 x 
10 4 mjy deg -2 to the 850 /im EBL if we exclude the A2390 
field because of the anomalously high number of red ob- 
jects, and 2.59 ± 0.19 x 10 4 mjy deg~ 2 if we include the 
A2390 field. These contributions are, respectively, 44 — 62 
percent and 59 — 83 percent of the 850 /im EBL (the ranges 
depend on the measurement of the EBL assumed) and are 
much larger than the 12 — 16 percent from galaxies con- 
taining bright AGN (Barger et al. 2001). ERO clustering 
is an important effect that must be accounted for if we are 
to obtain a better estimate. The median magnitude of the 
sources giving rise to the EBL is K' = 20 to K' — 20.3, 
depending on the cluster sample used. The 850 /im EBL 
contribution primarily arises in the subset of the popula- 
tion corresponding to dusty starbursts. 
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